Harmonic emission from cluster nanoplasmas subject to intense short laser pulses 
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Harmonic emission from cluster nanoplasmas subject to short intense infrared laser pulses is 
studied. In a previous publication [M. Kundu et al, Phys. Rev. A 76, 033201 (2007)] we reported 
particle-in-cell simulation results showing resonant enhancements of low-order harmonics when the 
Mie plasma frequency of the ionizing and expanding cluster resonates with the respective harmonic 
frequency. Simultaneously we found that high-order harmonics were barely present in the spectrum, 
even at high intensities. The current paper is focused on the analytical modeling of the process. We 
show that dynamical stochasticity owing to nonlinear resonance inhibits the emission of high order 
harmonics. 

PACS numbers: 36.40.Gk, 52.25.Os, 52.50.Jm 
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I. INTRODUCTION 

The study of rare gas and metal clusters interacting 
with intense infrared, optical, and ultraviolet laser pulses 
has emerged as a new promising research area in strong 
field physics (see [T], 0] for reviews). The generation of 
fast electrons and ions, the production of high charge 
states, the generation of x-rays and nuclear fusion (in the 
case of deuterium-enriched clusters) was observed using 
cluster targets in strong laser fields of intensities up to 
10 20 W/cm 2 . Hot and dense nonuniform, noncquilibrium 
and nonstationary plasmas produced under such condi- 
tions and confined on a femto- or picosecond time scale to 
nanometer sizes — so-called nanoplasmas — are new phys- 
ical objects with unusual properties. 

One of the most important features of nanoplasmas is 
the very efficient energy transfer from light to charged 
particles, which is much higher (per particle) than for 
an atomic gas of the same average density [3]. Al- 
though the particular mechanisms responsible for the 
laser energy deposition in clusters still remain debated 
i, S E H H ©, E9, HO, El, the pivotal role of collec- 
tive plasma dynamics (in particular the excitation of sur- 
face plasmons) and of nonlinear resonance (a definition is 
given below in Sec. IIIip in the energy transfer from laser 
light to the electrons in the nanoplasma and the subse- 
quent outer ionization was proved in experiments, simu- 
lations, and simple analytical models [rl M fol. [Tol. Il3l. [l4| . 

While the energy absorption from intense laser pulses 
by nanoplasmas has been widely studied both in experi- 
ments and theory, much less attention has been paid 
to harmonic emission from such systems. Naturally the 
question arises whether laser-driven clusters, being highly 
nonlinear systems, may be a source of high-order har- 
monics as efficient as atoms in the gas phase (or even 
more efficient?). There are at least two essentially differ- 
ent physical mechanisms which could be responsible for 
emission of harmonics from such a system. 

First, the mechanism based on recombination of the 



virtually ionized electron with its parent ion [l5j . well- 
studied in gaseous atomic targets, may also work in clus- 
ters where it can be modified by the fact that the atoms 
are closer to each other, so that the electron's motion be- 
tween ionization and recombination can be distorted by 
the field of the other ions. Moreover, the electron may re- 
combine with an ion different from its parent. This leads 
to modifications in the single-electron dynamics and in 
the phase matching conditions while the physical origin 
of harmonic generation (HG) remains the same as in a 
gas jet. Modifications of the atomic recombination mech- 
anism in clusters have been considered in [l6l [l7| . How- 
ever, in high intensity fields where each atom is loosing 
one or several electrons already during the leading edge of 
the laser pulse the recombination mechanism can hardly 
be efficient. 

Second, as dense electron plasmas is produced inside 
a cluster, its coherent motion may be an efficient source 
of radiation. In a collisionless plasma, as it is gener- 
ated by intense laser fields inside small clusters, individ- 
ual electron-electron and electron-ion collisions cannot 
destroy the coherency of the collective electron motion. 
This coherent, collective electron motion may cause HG if 
it is nonlinear. Recently, several experiments on high har- 
monic emission from plasma surfaces illuminated by short 
laser pulses of intensity 10 17 W/cm 2 and higher were re- 
ported (see, e.g., [H|,[l9j]). The physical picture which is 
behind such plasma harmonics appears to be more com- 
plicated and diverse than the recombination mechanism 
in atomic HG. 

In macroplasmas the magnetic component of the 
Lorentz force is a typical source of nonlinearity. In this 
case the nonlinearity parameter is v/c where v is the 
typical velocity related to collective oscillations of the 
electron plasma and c is the speed of light. As a conse- 
quence, generation of harmonics from dense macroplsmas 
may require relativistic intensities (l8| . As compared to 
macroplasmas, clusters introduce an extra source of non- 
linearity due to their small spacial size, namely X/Rq, 
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where Rq is the cluster radius and X is the amplitude 
with which the electron cloud oscillates under the ac- 
tion of the laser field [20j ]. Therefore, one could expect 
strongly nonlinear electron motion even in the nonrel- 
ativistic regime. It should be noted, however, that al- 
though individual collisions may not be important, there 
are other effects which can spoil the coherency required 
for efficient radiation. Most of these undesirable effects 
can be attributed to dynamical instabilities induced by 
the interaction of particles with the mean self-consistent 
field (in the presence of the laser field). Therefore, the ex- 
amination of plasma harmonics usually requires not only 
the analysis of the collective motion but that of individ- 
ual electron trajectories as well. 

Up to now only a few ex per iments on HG from clus- 
ters are known. In Refs. [2 lL I22I |23| HG from rare- 
gas clusters irradiated by infrared pulses of moderate 
(~ I0 13 — 10 14 W/cm 2 ) intensity was studied. It was 
shown that under such conditions harmonics can be gen- 
erated up to higher orders and with a higher saturation 
intensity than in a gas jet. In addition other interesting 
properties including different power laws in the intensity- 
dependent harmonic yield for particular harmonics for 
clusters and atoms were reported [21j . However, in 
Refs. [U, [12, HH the applied intensities were rather low to 
create a dense nanoplasma inside clusters and the main 
features of the recorded spectra (like a plateau followed 
by a cutoff) were found to be quite similar to the case of 
a gaseous target. This allows to attribute the observed 
effects to HG along the standard atomic recombination 
mechanism as described above. 

Recently the experimental observation of the third har- 
monic (TH) generation from argon clusters in a strong 
laser field was reported [24[ . The laser intensity was 
varied between 10 14 and 10 16 W/cm 2 and was thus suf- 
ficiently high for the creation of a dense nanoplasma. 
At such intensities HG along the atomic recombination 
mechanism relevant in gases is essentially suppressed be- 
cause of saturated single-electron ionization so that the 
observed TH signal can be fully attributed to the non- 
linearity of the laser-driven nanoplasma. A resonant en- 
hancement of the TH yield (when the Mie frequency of 
the expanding cluster approaches three times the laser 
frequency) has been measured using a pump-probe setup. 
The TH enhancement of the single-cluster response, as 
studied in theory before [2(| HH, [26| , is, however, masked 
in the experiment by phase matching effects whose op- 
timization at high average atomic density necessary to 
create clusters has been shown to be more intricate than 
for rare gas jets. The latter complication makes an ex- 
perimental study of nanoplasma radiation a difficult task 
while in computations it can be simplified by first exam- 
ining the single-cluster response and, second, analyzing 
propagation effects. 

In the recent paper (27[, we considered the radiation 
emitted by a single cluster exposed to a strong laser 
pulse. We computed harmonic spectra from argon clus- 
ters in short 800-nm pulses of intensity 2.5 x 10 14 — 7.5 x 



10 17 W/cm 2 using a 3D particle-in-cell (PIC) code (ap- 
plied before to the study of collisionless energy absorp- 
tion in laser-driven nanoplasmas [10]). The most intrigu- 
ing outcome of our study was the absence of high-order 
harmonics in the computed spectra, even at high inten- 
sities for all cluster sizes we considered. We attributed 
this effect to the above-mentioned dynamical instability 
in motion of individual electrons. 

In this paper we introduce two analytical models which 
describe the collective and single-electron dynamics of a 
laser-driven nanoplasma, respectively, and apply them 
for the explanation of both the numerical results [271 ] and 
the data [24[ . We show that the numerical results can be 
well described and understood within the rigid sphere 
model (RSM) but only for low harmonic orders 3,5,7 
while for higher harmonics the RSM yields qualitatively 
wrong predictions. Using a simple ID model we describe 
a stochastic, resonant single-particle electron dynamics 
which suppresses the emission of high harmonics. 

The paper is organized as follows. In Sec. [IT] we for- 
mulate the statement of the problem and describe the 
numerical method and the PIC results. In Sec. IIIII the 
spectra extracted from the PIC simulations are compared 
with the predictions of the RSM. In Sec. lIVI we introduce 
a model for the description of the single-electron dynam- 
ics and radiation and use it to explain the suppression of 
the high harmonic yield. The last section contains the 
conclusions. 



II. STATEMENT OF THE PROBLEM AND 
PREVIOUS RESULTS 

A cluster is converted into a dense electron plasma al- 
most promptly if the laser pulse is intense enough to ion- 
ize the cluster constituents. We refer to this process as 
inner ionization which should not be confused with outer 
ionization, when the nanoplasma electrons leave the clus- 
ter. Several competitive processes govern the evolution 
of this plasma during the interaction with the pulse and 
later, until the cluster becomes dissolved due to Coulomb 
explosion or hydrodynamic expansion. The electron den- 
sity increases because of further ionization of atoms and 
ions by the local electric field which may differ essentially 
from the applied laser field. With increasing plasma den- 
sity due to inner ionization, the oscillating electric field of 
an infrared laser is screened, so that its amplitude inside 
the cluster may be a few times or even an order of mag- 
nitude less than the amplitude of the incident wave (see 
more explanations on screening in Sec. IIVI below Eq.([8])). 
Simultaneously, as soon as a sizeable fraction of electrons 
have left the cluster, a quasistatic space charge field is 
built up which may become strong enough to induce fur- 
ther inner ionization ( "ionization ignition" , [2^, [23] ) . On 
the other hand, both outer ionization and the expansion 
of the cluster reduce the electron density. The net result 
of this competition is very sensitive to all parameters, 
including laser intensity, pulse duration, cluster size and 
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type of atoms. However, for the vast majority of param- 
eters a significant part of the electrons remains confined 
within the expanding ionic core. During this stage of 
the cluster evolution until the laser pulse is off the non- 
linear motion of the laser-driven nanoplasma may cause 
emission of laser harmonics. 

We have observed the above described scenario in PIC 
simulations of laser-cluster interaction, as reported in 
[ToL l27l ]. The dynamics and the radiation of nanoplas- 
mas was studied for Ar^r clusters (with the number of 
atoms N as 10 4 -10 5 and radii i?o ~ 6-10 nm), irradi- 
ated by linearly polarized, 8-cycle sin 2 -laser pulses with 
an electric field E;(t) = Eo sin (a;/i/2n) cos(u>it) and the 
wavelength A = 800 nm. Here u>i is the carrier frequency 
and n = 8 is the number of optical cycles in the pulse. 
The contribution of electrons remaining bound in atoms 
(ions) during the interaction was not taken into account. 
Results were reported in Ref. [27j where also more details 
about the numerical simulation may be found. Here we 
restrict ourselves to briefly summarize the main results: 

(i) The relative yields of low order harmonics depend 
on the laser intensity and the cluster size. 

(ii) The time-frequency analysis of harmonic spectra 
shows that low order harmonic enhancements occur 
when multiples of the laser frequency resonate with 
the transient Mie frequency. 

(iii) Even for the very high intensity 7.5 x 10 17 W/cm 2 
no distinguishable high harmonics (higher than the 
7th) appear in the spectra (see Fig. 1 in Ref. [27|). 

(iv) Only a part of the nanoplasma deeply bound inside 
the ion core contributes to (low order) harmonic 
generation. 



III. RESONANT ENHANCEMENT OF 
LOW-ORDER HARMONICS 

The enhancements of particular low-order harmonics 
predicted in [2(| HE] have been later studied numerically 
for small clusters using molecular dynamics simulations 
[26| and were finally observed in the experiment for the 
case of the third harmonic [H[ . It is well established now 
that the physical origin of enhancements is the resonance 
between the harmonic frequency and the Mie frequency 
of the expanding nanoplasma. The time-frequency anal- 
ysis of the radiation as calculated from the PIC results 
for the total acceleration further confirmed this state- 
ment [13] • Typical TF diagrams, showing which frequen- 
cies are emitted when, are reproduced in Fig. [T] In the 
same plots we show the scaled time-dependent effective 
Mie frequency UMie(t)/<-L>i- The standard definition of the 
Mie frequency WMio = \J 47rzno/3, where z is the average 
ion charge and n$ is the atom density, is only appropriate 
for the case of an almost homogeneous charge distribu- 
tion (in this section atomic units h = m = e = 1 are 



used). However, because of the ignition effect and the 
cluster expansion the average ion charge z depends on 
the ion position. As we found in [23], the low-order har- 
monics are emitted mainly by electrons confined within 
the ion core and with excursion amplitudes comparable 
to the initial cluster radius Rq or less. In the follow- 
ing we will refer to these electrons as the deeply bound 
electrons, and we define the effective Mie frequency as 
^Mie(i) = \j Qo{t)/Rf> wrt h Qo(t) the total ionic charge 
inside the sphere of radius R within which the cloud of 
deeply bound electrons oscillates. 

It is clearly seen from the TF diagrams of FigQ] that 
the enhancements of the third and fifth harmonic must 
be attributed to the respective high-order nonlinear res- 
onances between the laser frequency and the Mie fre- 
quency, u>Mie(t) = muj i (with m — 3, 5). Here we call the 
resonance nonlinear since it appears in a nonlinear sys- 
tem. This nonlinearity (i.e., anharmonicity in the effec- 
tive potential of electron-ion interaction) is a necessary 
condition for the emission of harmonics. The physical 
origin of this nonlinearity may be cither the electrons 
jutting out of the core during their motion, thus sensing 
the Coulomb tail [2(J, or the inhomogeneous charge dis- 
tribution within the ion core |25| . Note that the appear- 
ance of a significant signal at frequencies different from 
odd multiples of the laser frequency (namely, the second 
harmonic in Figflja) and the forth harmonic in FigOJc)) 
should be attributed to the excitation of eigenoscillations 
of the electron cloud with the time-dependent eigenfre- 
quency u>Mie(t)- These oscillations are of significance only 
in short pulses, as used in the PIC simulations. 

The results of the TF analysis in Ref. [27j not only ex- 
plain, at least qualitatively, the resonant enhancement of 
the TH observed in [24j but also confirm the idea formu- 
lated there that with a proper adjustment of parameters 
harmonics of higher order (5th, maybe 7th) could also be 
resonantly enhanced. 

It was shown in earlier studies [1(| HH that the reso- 
nant enhancements of low-order harmonics may be rea- 
sonably described using a very simple rigid sphere model 
(RSM). Here we show that, while for low-order harmonics 
the RSM works well, it fails even qualitatively to repro- 
duce the high-energy part of the spectrum. In a RSM 
@, [l(| H3] it is usually assumed that both ions and elec- 
trons form homogeneous, rigid spheres with sharp bound- 
aries. In this case the electron and ion charge density 
distributions are 

Pe(i)(r) = Tzn 6(R - r), (1) 

where no is the atom density in the cluster and 9{x) is 
the Heaviside step function. The signs =F correspond to 
electrons and ions, respectively. Within this model the 
restoring force F e j = uo^R^i depends upon the displace- 
ment X = i?oY of the electron cloud as [30] 

f(y) = — 9o(y), 
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troduced. The equation of motion reads 
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where E/((/3 = u;zt) is the electric field of the laser wave 
and 7 is the effective damping constant which can be 
estimated assuming a collisionless damping mechanism 
[H, [l2j (in calculations we use 7 = 0.1, in accordance 
with such estimates). 

Using the RSM significant insights into the absorption 
mechanisms in clusters and TH emission [2(| have 
been gained. However, an essential shortcoming of the 
RSM in the form specified by Eqs.([T|),(r5]) is the appear- 
ance of even powers of y in |f2| at y < 2. Physically it is 
the result of the discontinuous charge distributions (pj at 
r = Rq. As a consequence, in a perturbative treatment 
of harmonic generation the intensity of the s-th harmonic 
is proportional to E 2s ~ 2 instead of E 2s . For example, for 
the third harmonic one obtains with the RSM P 3uJ ~ Eq, 
which is obviously unphysical [20| . To improve the model 
it is sufficient to assume a smoothed charge distribution 
for the electrons while the positive charge density still 
may be described by (fT]). We use a Gaussian of charac- 
teristic width R e for the electron charge density, 



p e {r) = -zn exp (-r 2 / R 2 ) 



(4) 



and assume that the net charge density is zero at the 
cluster center. From (f4|) the total number of electrons is 
N e = ir 3 / 2 zn R 3 . The restriction N e < zN gives R e < 
(4/30i : ) 1 ^ 3 ^o ~ 0.91i?o- The equation of motion has 
now the same form (f3]) but with a modified restoring 
force: 



9o(y) -> gi(y) 



2 9 

ujfy 2 



(y 3 + l)erf [ofo + 1)] - (y 3 - l)eri[a(y - 1)] 



[l-2a 2 (l + j/ + y 2 )]- 



FIG. 1: (Color online) Time-frequency diagrams for an 
Ari7256 cluster at laser intensities 2.5 x 10 15 W/cm 2 (a), 
7.5 x 10 17 W/cm 2 (b). Panel (c) shows the TF diagram for an 
Ar92096 cluster at the intensity 7.5 x 10 17 W/cm 2 . The scaled 
time-dependent Mie frequency ajMio^)/^; is included in all 
plots (solid lines). 
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Here the dimensionless coordinate y and force f are in- 



,- a2 (y+V 2 [l-2a 2 (l-y + y 2 )] 



(5) 



where a — Ro/R e > 1.1 and erf (a;) = 
2/\/tt J q exp(—z 2 )dz is the error function. Con- 
trary to Eq. ([2|), a decomposition of the restoring force 
([5]) contains only odd powers of y while it has the same 
asymptotic Coulomb behavior for large displacements. 
The respective asymptotic expansions have the form 
[compare with (f2|)] 
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FIG. 2: Harmonic spectra (normalized to the signal at the fundamental frequency) calculated within the RSM for a cluster of 
radius Rq = 6.2nm, a — Ro/R e = 2, and max[oj c ff] ~ 3.2ui (a,d), max[w e ff] w 4Aui (b,e) and max[oj c ff] w 5.4a;; (c,f). The laser 
intensities are given in the panels in units of 10 16 W/cm 2 . 



with Ai — erf(a) — laj^fii exp(— a 2 ) > and A3 = 
— 4a 5 /5v / 7r exp(— a 2 ). It is clear that A\ — > 1 if a — ► 00, 
i.e., when the tail of the electron density distribution does 
not stick out of the ion core. In this case the eigenfrc- 
quency of small oscillations y^i^Mic is equal to o->Mic- 
For a finite a the spread of the electron cloud beyond the 
ion core reduces the eigenfrequency. 

Note that within the model we describe above the 
restoring force is nonlinear because a part of the electron 
cloud spreads out of the ion core. In this case the leading 
term in the nonlinear part of the force can be estimated 
as Uy ie XI 2 / RoR e - As a result, the third harmonic sig- 
nal scales as Rq with the cluster radius, in contrast to 
the standard Rq low expected for a volume bulk effect. 
This shows that the above discussed low-order harmonic 
emission from the nanoplasma is a surface effect which 
becomes relatively less pronounced with increasing of the 
cluster size. In small cold metal clusters subject to in- 
frared laser pulses of moderate intensity the Rq depen- 
dence was observed in [3l| though the experiment was 
not entirely conclusive (see also 32] for the respective 
theoretical analysis). 

Figure [2] shows harmonic spectra P(w) calculated from 
the above-described RSM. Within the model, the dy- 
namics and radiation of the electron cloud are governed 
mainly by the values of the two parameters lumic/^i and 
Eq/uj^Rq, determining the possibility of resonant en- 
hancements for particular harmonics and the number 
of excited harmonics in the spectrum. The eigenfre- 
quency of the electron cloud w e ff depends upon the ampli- 



tude of the oscillations (the maximum value max[a> e ff] = 
^/AiujMie corresponds to small harmonic oscillations) 
and, therefore, upon the laser intensity. As a result, at 
fixed cluster parameters the resonant enhancements ap- 
pear at certain values of the intensity. For the spectra 
of FigJ5][a,d) and (c,f) the parameters were chosen such 
that the maximum value of the eigenfrequency exceeds 
the integers 3 and 5, respectively, so that, with increas- 
ing intensity resonant enhancements of the third and the 
fifth harmonic appear. These enhancements are clearly 
seen from the comparison of these spectra with the ones 
calculated for the case shown in Fig(Ub,e) when the max- 
imum eigenfrequency max[w c g] ~ 4.4 is far from both res- 
onances. In panels (a-c) the intensities are chosen such 
that the sphere's oscillation amplitudes are almost the 
same in each case so that the relative difference in the 
strength of the third and the fifth harmonic (about three 
orders in magnitude) shows the efficiency of the resonant 
enhancement. For the higher intensities in panels (d-f) 
these enhancements are even more pronounced. 

A comparison between the RSM and the PIC results 
shows that for the low-order harmonics the RSM pro- 
vides a qualitatively correct description of the spectrum, 
including the effect of resonant enhancement. This gives 
an additional justification of the RSM for applications to 
laser-driven nanoplasmas. However, for relatively high 
harmonics (ninth or higher) the results predicted by the 
RSM appear to be qualitatively wrong. Indeed, the RSM 
predicts the appearance of higher harmonics with increas- 
ing laser intensity, as it is seen from Fig. f^f) where the 
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spectrum ends up at the 13-th harmonic. By a moder- 
ate variation of the parameters inherent to the model, 
namely u)y\\ e /uj and a, one may obtain even higher har- 
monics for the same intensities used for the spectra of 
Fig. [2] As expected for a classical system, no signature of 
a plateau emerges in the spectra. However, the PIC simu- 
lations reported in [53] do not show any harmonics above 
the seventh within the studied domain of parameters, 
including the very high intensity of 7.5 x 10 17 W/cm 2 . 
The same effect was found in |13| where no harmonics 
above the ninth have been observed for intensities up 
to 10 17 W/cm 2 . In Ref.[27j] we concluded from the in- 
spection of individual PIC electron's trajectories that a 
dynamical instability induced by the resonant interac- 
tion of electrons with the time-dependent self-consistent 
field is responsible for the suppression of high harmonics. 
Obviously, such a mechanism is beyond the RSM since 
the latter accounts for collective electron dynamics only. 
In the next section we introduce a single-particle model 
suited for the analysis of the motion and the radiation 
of individual electrons so that the instability responsi- 
ble for the suppression of high harmonics from cluster 
nanoplasmas under the above-described conditions can 
be studied. 



IV. SUPPRESSION OF HIGH HARMONICS 
VIA STOCHASTIC BEHAVIOR OF INDIVIDUAL 
ELECTRONS 

It is known from numerical studies (see examples in 

d m us Hi) that in a laser-driven cluster the elec- 
tron population separates into a dense core with a radius 
comparable to the initial cluster radius Rq and a rarefied 
halo with a typical size of several Rq. This subdivision 
is equivalent to a separation of quasifree electrons into 
deeply and weakly bound electrons, correspondingly. A 
decrease of the electron density due to outer ionization is 
compensated by the inner ionization of atoms, provided 
inner ionization is not yet depleted, so that the cycle- 
averaged electron density distributions both in the core 
and the halo evolve rather slowly in time. The core oscil- 
lates with relatively small deformations so that the RSM 
seems to be applicable to the deeply bound electron's dy- 
namics while electron trajectories in the halo are strongly 
disturbed by the laser field and cannot be captured by 
the RSM. 

In Ref. [27j the radiation of individual PIC electrons 
has been considered (see Fig. 3 there). It was shown 
that electrons radiate harmonics as long as they move 
inside the dense core. Being liberated from the core, 
electrons leave the cluster vicinity almost promptly, usu- 
ally within a laser period. This means, that the halo 
contains basically no permanent population, but consists 
almost entirely of electrons on their way out of the clus- 
ter. Hence the electron density in the halo determines 
the rate of outer ionization from the cluster. During the 
ejection, each electron emits an intense flash of radia- 



tion with an almost continuous spectrum that extends 
up to significantly higher frequencies than present in the 
net harmonic spectrum. In Ref. [13] we argued that these 
flashes add up incoherently in the total emission ampli- 
tude. Here we introduce a simple analytical ID model 
which helps to illuminate the physical origin of this in- 
coherency. Despite its simplicity our model is able to 
describe, at least qualitatively, all essential features seen 
in the simulations. 



A. Model 

Let us suppose frozen ions and the number of 
nanoplasma electrons fixed so that in the absence of the 
laser field each electron moves in the time-independent 
self-consistent potential U(x). The laser field excites os- 
cillations of the electron cloud which induce the ac part 
of the space charge field EsJ (t) . The net oscillating field 
inside the system can be written as 

£{t) =E l {i)+Eg\i) wfo/ttWwjt + a), (7) 

where f(t) is the time-dependent pulse envelope. If 
the Mie frequency notably exceeds the laser frequency, 
WMic > wi, as it is the case for infrared lasers, the laser 
field and the ac space-charge field E^J (t) essentially com- 
pensate each other. In this case the amplitude of the net 
oscillating field inside the cluster £ q is related to the am- 
plitude of the laser field E according 



LofE 



J Mic 



«£o- 



(8) 



This result is commonly referred to as screening of a low- 
frequency laser field inside clusters (see, e.g., the review 
0). There is no contradiction between this screening 
of the laser field for individual electrons and the fact 
that the whole electron cloud feels the unscreened field. 
Indeed, in the RSM there are two forces acting on the 
electron cloud: one due to the interaction with the ion 
core, another one due to the laser force, see Eq.®. In 
the single-electron description we should also take into 
account the interaction between the electron under con- 
sideration and all other electrons in the cloud. Within 
the model we assume that the electron cloud undergoes 
small, slightly nonlinear oscillations, so that this extra 
force is almost homogeneous inside the cluster and os- 
cillates in time with the frequency u>i. Within the RSM 
and under the condition uji <C uiMie we assume through- 
out the paper that the electron cloud displacement X(t) 
reads 



X(t) « - 



e£ 



m ( w Mi. 



LO 2 ) 



f(t) cos(w/< + a). 



Calculating the electric field induced inside the cluster 
due to this displacement and summing it up with the 
laser field ((7|) one obtains the estimate © for the ampli- 
tude of the net oscillating field. This type of screening 
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results from the coherent superposition of the applied 
and the self-consistent field and has nothing in common 
with damping of electromagnetic waves in macroplasmas. 
The latter occurrs on the spatial scale of the skin-depth, 
which is in general much bigger than the typical cluster 
size we consider. 

Taking a cluster consisting of N w 1.7 x 10 4 Ar atoms 
(Rq w 6.2nm) with the average ion charge ~z ~ 6 and 
the degree of outer ionization r] ~ 0.5, one can esti- 
mate fi^Mic ~ 6eV. In a Ti:Sa laser pulse of inten- 
sity 5 x 10 17 W/cm 2 (E w 3a. u.) the amplitude of 
the oscillating held inside the cluster is according to (JSJ 
So ~ 0.2a.u., i.e., more than one order of magnitude be- 
low the amplitude of the applied laser held. The qua- 
sistatic part of the space-charge field E^J^ — E sc — E^J 
which traps electrons within the ion core can also be es- 
timated for the assumed values of r] and z. Namely, the 
field near the cluster edge is Eii 1 ^ w r/Nze/RQ w 3a.u., 
i.e., more than one order of magnitude above the oscil- 
lating field amplitude. From this estimate we conclude 
that the oscillating field inside the cluster usually remains 
small compared to the quasistatic space-charge field. 

Within the model the electron's evolution is governed 
by the Hamiltonian 



H(p,x,t) = ^— + U(x)-e£(t)x 
2 m 



H (j>,x)-e£(t)x (9) 



V 



+ e£(t) = -eEgp (x) + e£(t), (10) 



and the corresponding Newton equation 

- - 

dx 

where m and e are the electron mass and charge. The 
well U (x) is created by the quasistatic part of the space 
charge. We model it by the function 



U(x) =t7 [l-l/Vl + (x/R Q n (11) 

where the values Rq and Uq are the cluster radius and 
the depth of the self-consistent well, respectively. Here 
we choose the energy minimum e = 0, so the e = Uq is 
the continuum threshold. For small excursion amplitudes 
this well is a nonlinear oscillator while for large excursions 
it has the desired Coulomb behavior. According to the 
estimates given above we assume the inequality 

e£n , ^ Uq 

« 1 - -^o = TT 
/to 

being satisfied, where Tq has the meaning of a character- 
istic quasistatic force trapping the electron. 



J~Q 



(12) 



B. Dynamics 

The dynamics of the unperturbed system with the 
Hamiltonian Hq is characterized by the energy depen- 
dence of the eigenfrequency 



n( e ) 



2tt 

W) 



T(e) 




(13) 
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FIG. 3: Dynamical characteristics of the potential the 
scaled eigenfrequency fi(e)/fi(0) (solid line) and the nonlin- 
earity parameter (|14|l (dashed line) vs the scaled energy e/t/o- 
Positions of the first and third-order resonances are shown by 
dotted lines. The value of the parameter j3 at the resonant 
energies is /3(ei) « 1.2 and w 4.5, respectively. 



where T(e) is the oscillation period, a(e), b(e) are the 
turning points, and e > is the total energy [34]. The 
energy-dependent parameter 



/3(e) 



(14) 



characterizes the nonlinearity of the unperturbed system 
and thus its potential capability to emit harmonics. Fig- 
ure [3] shows the energy dependence of the scaled eigen- 
frequency f2(e)/f2(0) and the parameter (fl"4)) . In cluster 
potentials the period T(e) increases with increasing en- 
ergy so that 



dn/de < 0. 



(15) 



The well pT]) has this property. 

The system (fl"0"| is not integrable. However, if the 
time-dependent force remains small in comparison to the 
first term at the right hand side of (fTU|) [i.e., under the 
condition (fT2|) ] the related dynamics can be qualitatively 
described using general methods of nonlinear mechanics 
[35l [H, [I?} ■ It is well established that the behavior of a 
nonlinear conservative system subject to a weak oscillat- 
ing field is essentially governed by the resonances between 
the energy-dependent frequency of the unperturbed os- 
cillations (fl~3|) and the frequency u>i of the external force. 
The positions e srn of these resonances are determined by 



sfl = mioi 



(16) 



with m and s integers. In a harmonic oscillator only 
the resonance with m = s = 1 can be excited. Reso- 
nances with m = 2,3,... are nonlinear with respect to 
the external field £(t) (in a quantum picture such reso- 
nances correspond to multiphoton absorption and emis- 
sion). Resonances with s = 2,3,... are related to the 
nonlinearity inherent to the unperturbed system Hq so 
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that for a highly nonlinear system they may be signifi- 
cant even in the weak perturbation regime. This situa- 
tion corresponds to the case we consider so that below 
we assume m = 1, disregarding the field-induced nonlin- 
earities. The reader should not be confused by the fact 
that we disregard here exactly the same nonlincarities 
which lead to the emission of harmonics we study in the 
rest of this work. Although these weak nonlinearities are 
responsible for the emission they play no role for the sta- 
bility analysis we consider in this subsection. Because of 
inversion symmetry of the well only odd resonances 
with s = 1, 3, 5... are important. Before the pulse an elec- 
tron moves along an unperturbed trajectory, specified by 
the energy eo and the initial phase $o- If the energy is 
far from any resonance level e s and the condition (|12D is 
satisfied the oscillating field in (JTOJ) induces a small per- 
turbation and, as a consequence, the trajectory is only 
slightly disturbed by the field. Under the condition (fT5j) 
the first-order resonance s = 1 is the lowest one on the 
energy scale so that the vast majority of such nonreso- 
nant trajectories lies well below the first resonant energy 
t\. This "perturbative" regime of interaction survives 
until the time-dependent energy approaches the vicinity 
of the resonance, either due to an increasing amplitude 
£o or because of a higher initial energy of the particle. 
Near the resonance the same small perturbation yields 
a strongly disturbed electron trajectory. A qualitative 
description of such trajectories and the determination of 
the boundaries that separate the "perturbative" and the 
"resonant" domains is possible using the methods de- 
scribed, e.g., in Refs. [HI, [H, [l?]]- We introduce new 
canonical variables, namely the "action" I and the "an- 
gle" fl, which are defined in the standard way [34l l35l. l3rj] 



Z7T 



p(e, x)dx, 



(17) 



Here p = ^/2m(e — U(x)) is the electron momentum in 
the unperturbed system and S(x) is both the position- 
dependent reduced action and the generator of the canon- 
ical transformation. The integration in (|17p is performed 
along the closed trajectory with the energy e. Strictly 
speaking, this canonical transformation is defined for the 
unperturbed Hamiltonian Hq only when closed, periodic 
trajectories exist. In this case one obtains immediately 
the new canonical variable $ = f2(e)t + #o- The vari- 
ables /, $ provide an effective zero-order approximation 
for the construction of a specific perturbation theory ca- 
pable of describing near- resonant dynamics [H, [H, H3] • 
Using the new variables and assuming that the particle 
energy is sufficiently close to a resonant value, say e%, we 
may omit rapidly oscillating terms and obtain an approx- 
imate, time-independent resonant Hamiltonian 

HrM) = H (r)-uiI--£x 1 (I)coB(il>), (18) 

so that the corresponding system is integrable. Here X\ 
is the first Fourier-component of an unperturbed trajec- 
tory, X(I, -i?) = J2k x k(I) cos(fctf) and ^ = i? - wrf. Now 




2 4 
Angle y (rad) 

FIG. 4: The phase space of the Hamiltonian (I19[l . A separa- 
trix (thick solid lines) subdivides the phase space into domains 
of finite (2) and infinite (1,3) motion. Two typical trajectories 
are shown by dashed lines. Dotted lines indicate the positions 
of the thresholds ±P m ax determined in (1201) . 



the trajectories of the system in the new phase space can 
be found, at least in the form of integrals. Further sim- 
plifications are possible if we take into account that for 
a weak perturbation the deviation of the action from its 
resonant value h = I(ei) is small, i.e., \I — I±\ -C I\. In- 
troducing the new canonical variable P = we obtain 
the nonlinear oscillator Hamiltonian (with the negative 
effective mass): 



M 



>0, 



B = 



(19) 



with = dQ/de\ e — ei . The electron behavior can thus be 
qualitatively described as nonlinear oscillations in (P, tp) 
space, known as phase oscillations [35l l36l [37| . In the 
new canonical variables which have the formal status of 
momentum (P) and coordinate (if>) the phase space of 
the Hamiltonian (p~9j) splits into domains of finite and in- 
finite motion (see FigQJ. Finite motion corresponds to a 
particle trapped by the resonance, while particles moving 
infinitely do not intersect with the resonance. The sepa- 
ratrix of (|19p is a boundary between these two domains. 
The motion near separatrices is unstable, so that even 
a small variation in the initial conditions may entirely 
change a trajectory. As a result, particles approaching a 
separatrix may penetrate from one domain to another or 
may be trapped by a resonance. The parameters charac- 
terizing the motion of trapped particles are the maximum 
deviation from the resonant action I\ and the frequency 
of small phase oscillations: 




9 



In energy space the positions of the separatrices of the 
first-order resonance are determined by 

4 = e(/i±P max ) = ei±Aei, Aei ~ U J^-~. (21) 

V -Ko P 

If the energy intersects a respective threshold so that 
| e — ex | < Aei, the electron becomes trapped by this res- 
onance domain and experiences phase oscillations with 
the frequency and amplitude both proportional to JJi. 
Due to the appearance of a new time scale given by the 
frequency of the phase oscillations (|20p the motion be- 
comes aperiodic and highly nonlinear. It should be em- 
phasized that in weak fields the perturbation parame- 
ter fi is far from resonance and ^JJi close to it so that 
the near-resonant motion appears to be much more per- 
turbed than the off-resonant one. 

The interaction with an isolated resonance cannot lead 
to ionization since the electron remains trapped for an, 
in principle, infinite time (38j. However, the higher-order 
resonances lying above may come into play. As soon as 
the separatrices of neighboring resonances intersect, the 
particle, captured by the first resonance, may jump to 
the third, etc. Because the volume of the accessible phase 
space is increasing with energy this inter-resonance mo- 
tion will have a predominant direction, namely towards 
higher energies. This leads to a fast liberation of the 
electron from the system, known as stochastic ionization 
@; HE H3 • For realistic parameters of laser-cluster inter- 
action this overlap of resonances is realized with almost 
100% probability so that cases where the particle remains 
trapped by a resonance are rare while almost prompt 
ionization occurs as soon as the first-order resonance is 
reached. 

We visualize the above-described scenario by solving 
Eq. ([lTj]) numerically. In the calculation we take u>i = 
Ro = 1.0 and Uq — 5.0 so that even the amplitude 
£q = 1.0 (/i = 0.2) still corresponds to the weak-field 
regime, as defined above. These parameters are not ar- 
bitrarily chosen. Indeed, a solution of Eq.fTU]) with the 
well (fTTj) and the field ([7|) with a slowly varying envelope 
depends on four dimensionless parameters, /i, T$ / mio 1 
eo/Uo and a, the two last of them defining the initial con- 
ditions. To recalculate all the parameters for a real sys- 
tem we should assume some certain values for the cluster 
radius and the laser field frequency, which then define all 
other parameters. For typical values, say Rq = 5nm and 
fkoi = 1.55eV, and for the dimensionless parameters of 
Fig. O one may check that the resulting strengths of the 
quasistatic and the oscillating parts of the self-consistent 
field indeed correspond to the estimates given below 
Eq.([8]). The positions of the most important first and 
third-order resonances are t\ m 0A8Uq and £3 « 0.77J7o, 
respectively. A particle with the initial energy eo < ei 
starts its motion at x — and (p = toit = —40 when 
the field £(t) is negligibly small. Then the field (with a 
Gaussian envelope) increases, and the electron propaga- 
tion under the action of the full force is calculated until 
Lp = +40. By choosing different phases a of the field we 



model different initial conditions for the particle at the 
fixed initial energy e - The results are summarized in 
FigH 

Figure [5^ corresponds to the "perturbative" regime of 
interaction. The initial energy is far enough from the 
first (lower) resonance, so that the trajectory in the en- 
ergy space does not intersect the respective lower sep- 
aratrix, or just touches it. As a result, the trajectory 
remains weakly disturbed, its shape is well described as 
a superposition of oscillations with the frequencies O(eo) 
and uji- By choosing different initial conditions we ob- 
tain trajectories simply shifted in time by the value of 
a. Figures EIb,c) correspond to the "resonant" regime of 
interaction where the initial energy is high enough or the 
field is strong enough to cause penetration of the particle 
into the vicinity of the first resonance. The five trajec- 
tories plotted in Figs[5jb,c) show that the near-resonant 
motion is very sensitive both to the initial conditions and 
to the field amplitude so that a particular trajectory ap- 
pears to be unpredictable. Usually the particle is emit- 
ted from the system while in rare cases it remains in a 
bound state after the pulse is off, being trapped by a 
resonance (see the trajectories in Fig. EIb,c)). From this 
observations we may conclude that at parameters typi- 
cal for intense laser-nanoplasma interactions the particle 
behavior in the "resonant" regime becomes stochastic, 
as it is expected to be according to the general theory 

[HUE If. 

In classical ionization, a particle has to overcome the 
potential barrier, i.e., its total energy must exceed the 
maximum value of the potential energy suppressed by the 
field at some time instant. Obviously, the trajectories of 
Fig. Etb,c) satisfy this condition, while both trajectories 
of Fig. do not. In Fig. [6] we show the trajectories 
in energy space evaluated for half the laser frequency, 
lu'i = lui/2 = 0.5, and all other parameters the same as in 
Fig. Eb- The plots show that with decreasing laser fre- 
quency the time-dependent energy gain from the field to 
the particle decreases down to a perturbative level, the 
total energy remains always below the maximum of the 
time-dependent potential energy, and no ionization or ex- 
citation occurs. Within the resonant picture we exploit 
here the qualitative difference between the trajectories of 
Figs. EIb,c) and [5] appears because with decreasing fre- 
quency the first resonant level is shifted up and no pen- 
etration into the resonant area between the separatrices 
takes place anymore. 

It is instructive to show the connection between the 
single-electron model specified by Eqs. (|7| . (|9|l . (fTOf . (fTT|) 
and a model which describes collective motion, as the 
RSM of Sec. HIT] does. The RSM deals with the electron 
cloud displacement X(i) whose Fourier-transform is di- 
rectly related to the spectrum. This value can also be 
calculated within a single-electron picture as 

X(i) = J de J da r(e ,a,i)F(e ,a), (22) 

where r(eo, a, t) is the individual trajectory with the ini- 
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FIG. 5: Upper panel: electron coordinate vs dimensionless time, x(tp — uiit), calculated from the numerical solution of (|10p 
for Uo = 5.0, u>i = 1.0, and Ro = 1.0, assuming the initial conditions x(— 40) = 0.0, < a < tt/2 (black solid and dashed 
curves) and a — n (gray curves). Other parameters are: eo = 0.2[/o, a — 0.0, £o = 0.5 (a), eo = O.SUo, a — 1.0, £o = 0.5 
(b) and eo = 0.2£/o, So = 0.8 (c) with a = 0.0 for a solid and a — 0.05 for a dashed curve. The electric field £(u>ot) is shown 
in the insert. Lower panel: scaled electron energy vs dimensionless time for the parameters of the upper panel, a-c. Thin 
horizontal lines show the position of the first-order (ei w 0.48(7o) and the third-order (53 ps 0.77C/o) resonances. Dotted lines 
show positions of the time-dependent separatrices (|21[) for the first resonance. A thick horizontal line shows the continuum 
threshold e = Uq. 
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FIG. 6: Same as the lower panel of Fig. [5p but for half the 
frequency, uj[ = 0.5 and the duration two times longer, so that 
the dimensionless time cp = uj[t varies within the same limits. 



tial energy eo and the initial condition a, and -F(eo, a) 
is the distribution function for electrons before the field 
is on. One should note that, contrary to the RSM, the 
spatial distribution of the electrons in the presence of 
the field is not known unless one calculates all individual 
trajectories. Calculating a trajectory r(eo,a,t) analyt- 
ically is possible only within perturbation theory with 



respect to the external field where one may easily derive 
© with the linear part of the restoring force only from 
([22]). The derivation of nonlinear corrections, although 
doable, requires very cumbersome algebra. The single- 
particle model is appropriate for a qualitative description 
of the stochastic resonant behavior but hardly applicable 
to the study of the slightly unharmonic motion of the 
deeply bound electrons. 



C. Radiation 

The analysis of the previous subsection gives a direct 
explanation of the radiation spectra extracted from the 
PIC results in Ref. |2jJ both for individual electrons and 
for whole clusters. Deeply bound electrons with energies 
below the first resonant separatrix e < t\ — Aei move 
along slightly perturbed regular trajectories. This causes 
HG with rapidly decreasing yield as a function of the har- 
monic order so that even the seventh harmonic is barely 
present in the corresponding spectrum of Fig. [7£i. An 
individual electron, while passing the resonance and be- 
ing trapped by it or leaving the cluster potential, emits 
radiation due to its strong acceleration, seen as a flash 
in the TF spectrograms of Fig. [7Jd,c. These spectro- 
grams should be compared with the ones extracted from 
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FIG. 7: TF diagrams corresponding to three (the black curve 
of Fig. [5^, and both curves of Fig. [5b) electron trajectories. 



our PIC results (see Fig. 3 in [27|). Exactly because of 
the stochastic nature of nonlinear resonance the elec- 
trons' trajectories are very sensitive to the initial condi- 
tions with which the nonlinear resonance is entered, as is 
clearly seen from Fig. [5fc where solid black, dashed black 
and gray trajectories correspond to the same initial en- 
ergy and the same field amplitude but different phases of 
the external field. As a result, flashes from different elec- 
trons are incoherent (the corresponding amplitudes have 
nearly random phases), and, being added up, vanish in 
the total dipole acceleration. 

This shows that exactly the same mechanism behind 
efficient energy absorption by and outer ionization from 
clusters, namely nonlinear resonance || [§, [l(| HH, re ~ 
stricts HG from them by breaking the coherent electron 
motion once it becomes strongly anharmonic. Only well- 
bound electrons trapped inside the ionic core with en- 
ergies far from the resonance contribute to the net, co- 



herent radiation of the cluster. A similar behavior was 
observed in classical ensemble simulations of atomic HG 

Instabilities induced by the resonant interaction grow 
in time (usually with an exponential rate). Thus the 
picture depends also on the pulse duration. With all 
other parameters fixed, for longer pulses, more and more 
trajectories from the vicinity of the separatrix experience 
a stochastic behavior. As a consequence, an increase of 
the pulse duration should lead, in general, to a further 
loss of coherency. 



V. CONCLUSIONS 

Although laser-irradiated cluster nanoplasmas emit 
low-order harmonics efficiently, no significant yield of 
high harmonics can be expected even for very high laser 
intensities. This is a consequence of dynamical stochas- 
ticity, inherent to nonlinear dynamical systems driven by 
weak, time-dependent forces such as the screened electric 
field inside a cluster. 

Increasing the laser intensity does not help much be- 
cause the self-consistent field trapping the electrons in- 
side the ion core increases too. As a consequence the 
electron population always splits into deeply bound elec- 
trons and a halo, the Mie-frequency and the screening 
increases, so that the physical picture remains almost in- 
sensitive to the intensity of the applied field. 

Another option we did not consider above is to use rel- 
atively long pulses where the first-order resonance with 
the Mie frequency can be reached because of the ion core 
expansion. If the Mie resonance is met the ac electric 
field inside the cluster is grossly enhanced and nearly all 
electron trajectories appear to be strongly disturbed. In 
this case an analysis of near-resonant stochastic behavior 
as given above is inappropriate. It seems that a direct 
numerical study is the only option under these condi- 
tions. An analysis based on classical Vlasov simulations 
was performed in Refs.[14]. 
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